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Abstract 

We investigate quantitatively the effects of strong electron-electron cou- 
pling onto the dynamics of lattice electrons. To that purpose the self- 
consistent version of the bubble-chain approximation at zero temperature 
and half filling of the Anderson (Hubbard) model is used. Special attention 
is paid to a critical region of an electron-hole correlation function shaping 
the transition from weak to strong interaction. We find an analytic solu- 
tion with Fermi-liquid properties on the weak-coupling side of the critical 
region around the two-particle pole. It is shown that Fermi-liquid theory 
does not lead to a consistent behavior of the self-consistent solution on the 
strong-coupling side of the critical region. 



PACS Numbers: 71.27+a, 71.28+d 
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The Anderson and Hubbard models provide a microscopic description of the effects 
of electron-electron correlations onto the dynamics of a lattice electron gas. Especially 
recently the single-impurity Anderson model (SIAM) has newly attracted attention of 
theorists because of its role in the exact description of the Hubbard model in d = oo |1]J2|]. 
However, the two models behave at strong coupling in qualitatively different manner. 
At half filling and zero temperature, we expect a Kondo-like behavior, i.e. a narrow 
resonance at the Fermi level, in the SIAM, while the Hubbard model in d = oo, when the 
antiferromagnetic LRO is suppressed, turns insulating. To understand the differences in 
behavior of these two related models, it is necessary to have an approximation reliable at 
intermediate and strong coupling for both Anderson as well as for the Hubbard model. 

Although we know much about the SIAM from the Bethe-ansatz solution [|3| this 
method has yet proved inefficient in the Hubbard model in d = oo in spite of an exact 
transformation of the d = oo Hubbard model onto a SIAM with a self-consistent condition. 
The only technique equally well applicable to the SIAM and to the lattice models is many- 
body perturbation theory summed via Feynman diagrams. 

We know from earlier studies on the SIAM |3]|| that only self-consistent (renormalized) 
sums of diagrams can provide reliable approximations at intermediate and strong coupling. 
Otherwise we cannot evade an unphysical RPA pole in a two-particle Green function |J . 

Renormalized sums of Feynman diagrams for the Hubbard model in d = oo at weak 
coupling were studied recently It was shown || that self-consistent, renormalized 

sums of the RPA-type can be used at any temperature and in principle at weak as well 
as at strong coupling. However, it is not straightforward to extrapolate such theories 
consistently to the strong-coupling regime. There is no analytic solution to these advanced 
renormalized sums and numerical solutions break down before the strong-coupling limit is 
reached. The numerical troubles arise when we are approaching the RPA pole (singularity) 
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in an electron-hole correlation function. There is then no effective way to make the 
iterations converge in the strong-coupling regime. It is then crucial to decide from analytic 
estimates how the pole in the two-particle function is approached by the full, numerically 
unreachable solution. 

The transition from the weak to the strong coupling regime can hence be investigated 
only analytically using some assumptions. First such an analytic study in the SIAM was 
done in ref. 0, where a low-frequency approximation was used to estimate the behavior 
of Suhl's renormalized RPA in the critical region of the two-particle pole. The aim of this 
paper is to reinvestigate the transition region between weak and strong coupling regimes in 
the SIAM and the Hubbard model in d = oo. We extend Hamann's approach from ref. || 
based on a Fermi-liquid, low-frequency expansion for electron-hole bubbles and show how 
dominant contributions to the self-consistent solution can analytically be estimated when 
the two-particle pole is being approached. As an example we use the self-consistent version 
of the bubble-chain (shielded interaction) approximation |7||| applied to the SIAM at half 
filling and zero temperature. This approximation, in contrast to earlier theories @|||9|], 
represents a thermodynamically consistent and conserving theory || applicable to the 
SIAM as well as to the Hubbard model. As a result we obtain a set of algebraic equations 
the solution of which forms a Fermi-liquid at weak coupling, but leads to inconsistent 
results at intermediate and strong couplings. This inconsistency is explained by the 
failure of the Fermi-liquid, low-frequency ansatz to capture all the relevant features of the 
full solution at intermediate and strong coupling. It is necessary to take into account also 
incoherent, non- Fermi-liquid contributions to the two-particle Green function to reach a 
consistent behavior of diagrammatic approximations at intermediate and strong coupling. 

The bubble-chain approximation for the self-energy can generally be written in the 
spin-polarized version as ]7||| 
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where X cr (iz/ m ) is a contribution due to the electron-hole bubble, 

^ oo 

^W = ^ XI G*( iv m + iu n )Ga{iUn) ■ (2) 
P n=—oo 

The electron propagator G a (z) is defined for the SIAM as G(z) = [z + p a — V 2 T a (z) — 
S^z)] -1 , where \i a = p. + ah — e/ is the effective chemical potential and T a (z) is the local 
element of the Green function of the conduction electrons. For the Hubbard model in 
d = oo we have G a (z) = J dep(e)[z + p — ^ a ( z ) — e ] \ where p(e) is the density of states 
(DOS). 

We can analytically continue the sums over the Matsubara frequencies (u n = (2n + 

l)7r/3 _1 ,z/ m = 2mirf3~ 1 ) to the real frequencies and after some manipulations we obtain 

the following representations at (3 = oo, n — 1 and h = 

o 

ReX(u + ) = J dxp(x)Re[G(x + uj+) + G(x - lu + )], (3a) 

— oo 

M 

lmX(uj + ) = —izsgmo J dxp(x)p(x — \uj\) , (3b) 



where p{uj) = — -lmG(u) + ), u + = u + i0 + . The self-energy can then be represented as 

U 2 } 

ReE(c<j + ) = — YJ^ X {p( x ^ e — UJ +) — C{x + co + )} 

— oo 

+ ilmC(a; + )Re [G(x - uu+) - G(x + , (4a) 

M 

Im£(u; + ) = U 2 J dxp{x - |u;|)ImC(x + ) . (4b) 
o 

The two-particle correlation function C{z) := X(z)/(1 — U 2 X(z) 2 ). Equations 
represent a set of nonlinear integral equations for Re£(u; + ) and Im£(u; + ). These equa- 
tions can be solved numerically by iterations at weak coupling but the iteration 



procedure breaks down as C(0) — > oo with increasing U. Since X(0) < the quantity 
1 + UX(0) approaches zero at intermediate coupling. The dominant contributions to £(u;) 
then come from a vicinity of the Fermi energy (u ~ 0) where the two-particle correlation 
function C(u>) is sharply peaked. We now use the Fermi-liquid assumption that only the 
low-frequency behavior around the Fermi energy is decisive for the physics of interact- 
ing electrons around the two-particle pole and replace the denominator of C(u) with a 
quadratic polynomial 

1 + UX(u) w UX"(0) [A 2 + u 2 - mau] (5) 

where X(Q)" := f _ oo dxp(x)ReG"(x), a := v 2 /X"(0), A 2 := (1 + UX(0))/UX"(0). Here 
v is the DOS of the unpertubed Green function at the Fermi energy. The parameter A is 
an energy scale measuring dominant fluctuations in the critical region of the two-particle 
pole. Note that Hamann used in || the same idea of a low-frequency expansion at the 
denominator of a two-particle function, but expanded X(u) only to linear power in u. 
This difference leads to drastic changes in the critical behavior of the solution. It is also 
necessary to realize that (g) is valid only if Fermi-liquid theory holds without restrictions, 
i.e. there are no other relevant energies except for the Fermi one. 
Inserting fl5|) in (|j) we obtain in leading order of the limit A — > 

ReS(cj + ) = ReG(u) J(oo) — p(u)sgnujK(uj) 

aReG(u) °r x p(cu)sgnu; A 2 + x 2 

2X"{0) J X {A 2 + x 2 ) 2 + 7c 2 a 2 x 2 ~ 2X"{0) J X {A 2 + x 2 ) 2 + n 2 a 2 x 2 ' ( & > 

nap(u) '/•' x 

ImS w + ) = — piu) J(uj) = — . . / ax— — ^— n - - . (6b) 

v +; HK ' y ' 2X"(0) J (A 2 + x 2 ) 2 + n 2 a 2 x 2 v ' 

We see that integral equations (|) turned algebraic, where only two positive parameters 
A and X"(0) are expressed as integrals over the products of the full Green function G(uj). 



The parameter X"(0) is proportional to the effective mass (—£'(())) of quasiparticles 
from Fermi-liquid theory and can be assumed as an effective mass of electron-hole pairs. 
The energy A is a new relevant scale for the two-particle scattering. Although (|6|) is 
strictly valid only in the limit A — > 0, we can extrapolate it also to the weak coupling, 
U — > 0, A — > oo. Such a theory then fulfills Fermi-liquid assumtions, i. e. ImE(w) ~ —uj 2 
and Re£(u;) ~ — uj as to — > 0, and ImE(w) ~ — 1/uj 2 and Re£(co>) ~ 1/u as uj — > oo. 
Approximation @ represents the simplest Fermi-liquid theory with frequency dependent 
self-energy determined essentially from algebraic equations. It may serve as an alternative 
to the recently proposed approximations trying to clarify the way Fermi liquid breaks down 
at strong coupling of the Hubbard model [|I|,p!0|-[T3|] . 

The integrals in (^]) can be performed explicitly. To simplify the studied equations we 
confine our analysis only to the SIAM. If we use the standard approximation T(u + ) = —iT 
we can resolve (0) analytically in the limit A — > 0. We find an explicit solution 



ImTi(uj) 



v±r 2 v 2 v 

1 — 2 H > ( 7a ) 

I \-\- ' J(oj)u — sgnu;i^(ct;)lmS(co' + ) ^ 2 



\ \J(uj)V 2 T - [J(u) + J(oo)]lmE(o; H 



ReE(o;) = -ImS(o; + ) — — - . 7b 

J(uj)V z T — [J{uj) + J(oo)JImL(c<j + ) 

To close the approximation we complete these equations with definitions of the parameters 



X"(0) 



dxp(x)ReG"(x) sa z/ReG"(0) , 
' 1 + 2U J dxp(x)ReG(x) 



(8a) 



b) 



C/Jf"(0) 

The set of equations (0) and @ can be solved numerically. Contrary to Hamann's result 
we reach a critical interaction U c ~ 3.7 /u at which A = and a pole in the electron-hole 
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Green function appears at the Fermi energy. The pole leads, however, to a quite different 
behavior of the self-consistent solution than in the non-self-consistent RPA. The reason 
for this behavior deviating qualitatively from that found in RPA on the one side and by 
Hamann on the other side lies in a breakdown of Fermi-liquid theory around U c . The 
expansion coefficient X"(0), neglected by Hamann, diverges at the critical point. Namely 
AX"(0) — > a c > and the effective mass of quasiparticles becomes infinite. Hence the 
low-frequency approximation @ indicates a sharp transition between weak and strong 
coupling regimes. The existence of a sharp transition is, however, incompatible with a 
strong-coupling solution, unless the critical point represents a metal-insulator transition, 
i.e. the DOS at the Fermi energy vanishes with A — > 0. Namely, it is easy to show that 
the solution at and above U c does not possess Fermi-liquid properties and 

Re£((T) = -ReS(0 + ) > (9a) 

and when |ReS(0)| < w, where w is a half bandwidth, then 

ImS(iO + ) < . (9b) 

The analyticity assumption for the expansion (|5|) hence does not hold any longer. 
Analyzing the equations at U > U c we find that (|5|) must be replaced by 

1 + UX(u) « UX'(0) [A + \u\ - iitau} (10) 

reflecting a nonanalyticity of the particle-hole bubble at low frequencies at strong coupling. 
The weak-coupling and strong-coupling ansatzes (|5|) and (|IUD, respectively, are evidently 
incompatible and do not allow a continuous matching. There is no critical point from 
the strong-coupling side as in Hamann's analysis. We must hence conclude that the 
above analysis based on @ is incomplete and does not lead to the actual strong-coupling 
asymptotics of self-consistent diagrammatic approximations. 
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Since the real part of the self-energy experiences a jump and the imaginary part acquires 
a nonzero value at the critical point, expansion (|^) around uj = in this form becomes 
meaningless. The frequency interval within which Fermi-liquid theory holds has shrunk 
to zero at U c . To assess the solution in the critical region, U /* U c , more precisely 
we have to take into account two contributions to the function X(u + ). The first one 
is that from (|5]) and holds now only for \u\ < O. Energy determines the interval 
within which Fermi-liquid theory holds. It can be defined as a frequency where the real 
(imaginary) part of the self-energy reaches its (first) extremum. It is essential that — > 
as A — - > 0. Function X(u+) outside the interval [—0,0] must be newly approximated 
and one has to expand around the points ±0. Such an expansion has its leading terms 
of type flltf). In the limit A — > the non- Fermi-liquid contributions from the additive 
expansion around — > more and more take over the control of the critical behavior 
and preclude the critical point to be reached. It means that the mechanism how the 
RPA pole is circumvented in self-consistent approximations with two-particle bubbles is 
much more complicated than anticipated in Hamann's low-energy analysis within Fermi- 
liquid theory. Only the full treatment of self-consistent diagrammatic approximations 
with the non-Fermi-liquid contributions fllCf ) can produce the genuine strong-coupling 
asymptotics. It is not yet clear whether the expected strong-coupling Kondo asymptotics 
in the SIAM will be reproduced correctly in this way. The result depends namely on a 
detailed balance of the Fermi- and non-Fermi-liquid contributions in the critical region 
and on the dispersion relation of the underlying lattice. The complete analysis of the 
bubble-chain approximation in the critical region A — > will be presented elsewhere. 

To conclude, we demonstrated that self-consistent diagrammatic approximations at 
intermediate and strong coupling, where a pole in a particle-hole correlation function is 
approached, show interesting behavior going beyond Fermi-liquid theory. We analyzed the 
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self-consistent version of the bubble-chain approximation for the Anderson and Hubbard 
models and proposed a simple analytic solution with Fermi-liquid properties at weak 
coupling. However, to obtain the asymptotics at strong coupling, it is necessary to restrict 
the validity of Fermi liquid theory to an interval around the Fermi energy, u G [— 0, 0], 
vanishing when the critical point is reached. The contributions obtained from Fermi 
liquid theory alone are insufficient to suppress the two-particle singularity and to reflect 
the strong-coupling behavior. Differences in solutions of @ and that of Hamann from ref. 
H show how a delicate problem it is to find the pertinent asymptotics of the Anderson 
and Hubbard models in the critical region of a two-particle pole crucial for the transition 
from the weak to the strong coupling. Neither of the above analytic solutions can yet be 
seen as exact at strong coupling within the chosen self-consistent approximations. 

The work was supported in part by the grant No. 202/95/0008 of the Grant Agency 
of the Czech Republic. 
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